On the theory of dielectric relaxation 
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Dielectric relaxation has been investigated within the framework of a modified mean field theory, in 
which two dynamic processes are believed to be involved in dielectric relaxation; the first corresponds 
to relaxation of crystalline lattice and the second represents slow relaxation of a partially ordered 
phase of high-temperature structure precursors generated by thermal fluctuations in solid dielectrics. 
It later becomes clear that it is the interaction between these two processes that results in a universal 
dielectric relaxation law. 

PACS numbers: 77.22.Gm,77.84.-s 



In 1913, Debye investigated the anomalous dispersion 
phenomenon, in which the corresponding index of refrac- 
tion falls with frequency, of a group of dipolar molecules 
[l| ; he treated a dipolar molecule of the group as a sphere 
immersed in a viscous fluid and then, under the assump- 
tion that the only electric field acting on the molecule 
is the external field, used Einstein's theory of Brownian 
motion to tackle the collisions between the molecule 
and its neighboring molecules in the liquid; eventually he 
formulated the following equation [l[, 
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where e is the complex permittivity of the dipolar 
molecules and u is the angular frequency of the external 
field; Eqo and £ s represent the permittivity at the high 
frequency limit and the static permittivity of the dipo- 
lar molecules, respectively; r c is the characteristic relax- 
ation time of the dipolar molecules. The above equation 
is often called the Debye relaxation law, which represents 
the dielectric relaxation response of the above-mentioned 
dipolar molecules to an alternating external electric field. 

The dielectric relaxation phenomena had been exten- 
sively investigated long before the Debye law was pro- 
posed. Perhaps it is the first relaxation law that was 
derived soundly based on statistical mechanics, thus 
the Debye law is often regarded as the starting point 
for studying relaxation responses of dielectrics. Un- 
fortunately, numerous experimental studies have shown 
clearly that the relaxation behavior of a wide range 
of dielectric materials deviates strongly from the De- 
bye law. Over the past 70 years, many empirical re- 
laxation laws, as the variants of the Debye law, have 
been developed. Among the most important are the 
Cole-Cole equation (1941-1942) 0,[5j], the Cole-Davidson 
equation (1950-1951)1, 0, the Havriliak-Ncgami equa- 
tion (1966-1967) d, the Kohlrausch- Williams- Watts 
function (the Fourier transform of the stretched expo- 
nential function or the Kohlrausch function [l(|) (1970) 
[Tlj . etc. In practice, these empirical relationships work 
very well for certain materials, but not for others. In 
1970s, Jonscher and his co-workers analyzed dielectric 
response data of around 400 hundred dielectric materi- 



als and then he suggested that there exists a universal 
law of the dielectric response [l2l - fl6j . Largely due to 
Jonscher's pioneering work, people have noticed an in- 
teresting fact that the relaxation behavior is very similar 
despite different materials used and experimental tech- 
niques employed and most relaxation data can be in- 
terpreted by two types of experimental fitting function 
[171 ]: the Kohlrausch function or the stretched exponen- 
tial function (1854) [l(| (historically, the citations related 
to the Kohlrausch function might be confused; the corre- 
sponding clarification has been given in Ref. [l8j), which 
is written below 



f(x) = cxp - 
or the Jonscher function [lj| 
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where r, a, and b are parameters for a given material. 
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FIG. 1. Diagrammatic representation of normal and distorted 
crystal lattice; (a) normal crystal lattice without net dipole 
moment; (b) and (c) distorted crystal lattice with net dipole 
moment represented by arrows. 

Surprisingly, the above-mentioned empirical laws or 
functions are all fractional power-law relationships. If 
a universal dielectric relaxation law does exist, then, at 
least, those laws or functions could be derived from the 
same fundamental equation. In this paper, we will at- 
tempt to find this fundamental equation and explore the 
physical mechanism underlying dielectric relaxation phe- 
nomena. 

As usual, we start with the Debye law. It is generally 
believed that neglecting the mutual interactions between 
dipolar molecules and their neighboring molecules is the 
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reason that makes the Debye law deviate from exper- 
imental observations. The Debye law, however, is the 
model based on statistical mechanics, which is closely re- 
lated to the Einstein theory of the brownian motion and 
the Smoluchowski equation [19j : like the collisions be- 
tween dipolar molecules and their neighboring molecules, 
the interactions have been implicitly considered in De- 
bye's model and the corresponding induced dipole mo- 
ments have been assumed to be stochastic in the time do- 
main and randomly distributed in the spatial domain so 
that the net influence of the interactions is zero. There- 
fore, the problem associated with the Debye law is not 
whether the interactions are neglected or not but what 
kinds of dipole moments can be induced by the inter- 
actions and cannot be averaged out. In other words, we 
have to take a close look at what kinds of dipole moments 
induced by the mutual interactions are neither purely 
stochastic in the time domain nor randomly distributed 
in the spatial domain and consider how they would be- 
have and then alter dielectric relaxation processes during 
the perturbation of external electric fields. Obviously, if 
we try to solve the problem, we must consider the in- 
volved many-body problem. In his recent voluminous 
work, Ngai has repeatedly emphasized that the universal 
dielectric relaxation law, if it exists, must be built on the 
basis of the many-body problem [20|]. Inspired by Ngai's 
profound insights on relaxation in complex systems, we 
will attempt to use a modified mean field theory (MMFT) 
to study the dielectric relaxation behavior. 

Low temperature phase ■ Nematic phase ■ Isotropic phase 



FIG. 2. Schematic representation of the cooperative behav- 
ior spectrum of HTSPs. Below temperature Tn and above 
temperature Tni, there is no cooperative movement in HT- 
SPs; between Tn and Tni, HTSPs can behave cooperatively 
to form a nematic phase. 

We now need to address how certain dipole moments 
could be induced and then explain what kinds of micro- 
scopic structures those dipole moments could form un- 
der an external electric field. For a crystalline material, 
it has a purely ordered structure at absolute zero and a 
completely disordered structure near its melting point; 
at temperatures between those two extremes, the mate- 
rial would possess both ordered and disordered structures 
largely due to thermal fluctuations (we here do not con- 
sider the defects and microscopic disordered structures 
generated during manufacturing processes in dielectric 
materials). Let's consider a crystalline dielectric lattice 
diagrammatically shown in Fig. [1]. Even at temper- 
atures far below the melting point, there always exists 
the probability that certain atoms in the normal crystal 
lattice shown in Fig. [la] could gain extra kinetic en- 
ergy to move quasi-permanently away from their original 
equilibrium positions and distort crystal lattice as shown 



in Figs, [lb] and [lc] due to thermal fluctuations. This 
kind of distortion can introduce the local non-uniform 
deformation, which breaks the inversion symmetry and 
induce dipole moments even in non-polar dielectric mate- 
rials [2l[; the induced dipole moments are schematically 
shown in Figs, [lb] and [lc]. The local disordered po- 
lar structures corresponding to such distorted crystal lat- 
tice are defined as high-temperature structure precursors 
(HTSPs). Since HTSPs possess higher potential energy, 
they are often metastable. Recent studies have shown 
that those polar HTSPs can form a partially ordered ne- 
matic phase over a wide temperature range from Tn to 
Tni due to the competition between energy and entropy 
in dielectric materials [22| . If we use v to represent the 
local preferred direction and approximately treat HTSPs 
as molecules, then their orientational order parameter, 
S op , can be written as (23l. |24| 

Sop = 1(3 (v\n) 2 1) = i((3cos 2 ^ - l)>, (2) 

where ( } represents the average; v l is defined as the 
given direction of the disordered structure located at the 
position i; n is usually called the director that represents 
a particular direction; 9 l is defined as the angle between 
v l and n at the position i. If < S op < 1, we can say 
that HTSPs cooperatively form a nematic phase. S op = 1 
corresponds to an ideal case, in which all HTSPs are 
perfectly aligned. If temperature continues to rise, at a 
certain point T = Tni, the thermal energy will be large 
enough to disturb HTSPs, which makes HTSPs randomly 
oriented and then renders S op = 0. Therefore, when T > 
Tni, all HTSPs behave like a normal liquid, which has 
an isotropic phase. The cooperative behavior spectrum 
of HTSPs is shown schematically in Fig. [2]. 
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FIG. 3. Diagrammatic sketch of the loading signal; the char- 
acteristic relaxation occurs only once within one T p , which is 
the period of the loading signal. 

Within the framework of MMFT, the Gibbs free energy 
of a unit volume crystalline dielectric specimen contain- 
ing polar HTSPs can be written as 

G = U-TS-(l + k)ED = U-TS-(l + k)^rP 2 , (3) 

or 

G = U-TS-(l-k)ED = U-TS-(l-k)^rP 2 , (4) 

where P and D represent the electric polarization and 
the electric displacement, respectively; Sq and e r repre- 
sent the electric permittivity of free space and the relative 
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permittivity, respectively; x is defined as electric suscep- 
tibility; E is the applied electric field; k here is defined 
as 



k S op {T - Tjv) 



T 



N 



(Tjv <T< T m ), 



(5) 



where both k and fco arc dimensionless coefficients. The 
physical meaning of k can be explained as following - it 
actually represents, in the statistical sense, the fraction 
of the potential energy, which is generated by the coop- 
erative movement of polar HTSPs. The term (1 + k) in 
Eq. ([3]) demonstrates that the crystalline and the ne- 
matic phases of the test specimen respond to an applied 
electric field, which often corresponds to the case that the 
field is too weak to disturb the nematic phase. However, 
the term (1 — k) in Eq. (j4j) represents the situation that 
is often encountered in practice; the negative sign here is 
due to Le Chatelier's principle, which means that, if the 
applied field is large enough to be able to disturb the ne- 
matic phase, the nematic phase as the ensemble of polar 
HTSPs would then undergo a specific structural change, 
in which those polar HTSPs would tend to huddle to- 
gether and grow up (this means that k will gradually 
increase when the external field is continually applied), 
to counteract any imposed polarization by the filed |25j . 
Clearly, in either cases, the evolution of k during the 
perturbation of the applied field could alter the linear 
polarization-electric field relationship or the dielectric re- 
laxation behavior of the test specimen. In order to deter- 
mine how the relaxation behavior is influenced by k 7 we 
need to exploit a modified Landau-Khalatnikov equation, 
which is given below [26| 
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dt 



dP' 



(6) 



where 7 is a kinetic coefficient and is considered to be in- 
dependent of temperature 0, H3| ; t represents the time 
variable. We will consider the following two cases sepa- 
rately. 

Case I - the applied electric field is weak and the term 
(1 + k) is given in the Gibbs free energy. Substitute G 
given in Eq. Q into Eq. ©, we get the following result 



at oP eox 



Solve this equation, we have, 
2(1 + k)e r 



P = P s exp 



= P s exp 


t~ 
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(7) 



(8) 



where P s is defined as the induced electric polarization 

2 

under static electric fields and r = 2 (\+k)s • ^ wc ^ ur ~ 

2 

ther assume that t c = 7 ^° Xe represents the characteristic 
relaxation time of the test specimen, then the relation- 
ship between r and r c can be written as r = or 
t = (I — k')r c , here k' = k/(l + k). Thus, t can be 



defined as the modified characteristic relaxation time of 
the test specimen. Clearly, the above equation can be 
written in the following two forms. 



or 



X = X s cxp 



s exp 



(9) 



(10) 



where both Xs and e s represent the static susceptibility 
and the static permittivity, respectively. 

The above equations, Eqs. ©, ©, and (JTDJ), only rep- 
resent a single modified characteristic relaxation after an 
instantaneous external perturbation. In order to investi- 
gate the dielectric relaxation behavior of the test speci- 
men, we need to determine the relationship between the 
characteristic relaxation and the applied field. For this 
purpose, we consider the applied electric field that has 
the waveform that is shown schematically in Fig. [3]. If 
defining T p as the period of the field, we can clearly see 
that it is a one-direction signal, in which there is no re- 
versal and, if T p 3> t, the characteristic relaxation will 
fully occur once at -f- within one period. The dielectric 
relaxation of the test specimen under an applied electric 
field reduce to a continuous and progressive relaxation 
process consisting of a series of the modified characteris- 
tic relaxations. Obviously, there is a high frequency limit 
imposed by the test specimen itself. If the frequency of 
the applied field reaches this limit or approaches r c , the 
specimen cannot fully respond to the field and its relax- 
ation will start to saturate. Without loss of generality, 
we define the value of e at this limit as £00. Then Eq. 
(fTUll can be modified as 
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We then take the Fourier transform of the right hand side 
of Eq. (fTTj). The result is given below, 
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cxp( )u(t) 
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1 + iuiT 



(lu = 2tt/), (12) 



here / is the frequency of the applied electric field. How- 
ever, t in Eq. (|11[) actually varies from to r not from 
to 00. We can imagine that we squeeze the waveform 
of the applied signal from (0, 00) to (0, r) in the time 
domain and, let's say, if r = 0, then the squeezed wave- 
form should behave like the Dirac delta function. We 
now modify Eq. (fT2"|) as follows 
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(13) 



where we can see that, when r — >• 0, (s(oj) — £qo)/(£s ~ 
£00) = 1; which is exactly the Fourier transform result 
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of the Dirac delta function. The above equation can be 
further expressed as 



e{u) 



1 + - k')r c 



(14) 



If k = 0, then k' = so that this equation will reduce 
to Eq. (fTJ), which is the mathematical expression of the 
Dcbye relaxation law 

If we assume that k' <C 1, then iuj(l — k')r c can be 
simplified as — k')r c al - 1 + k' + (1 — k')iu!T c — 

I + (1 — k')[iuT c — 1]; here wr c < 1 and, therefore, \iu)T c — 

I I < 1 . By taking advantage of Eq. (|F31|) in Appendix, 
we have iuj(l — k')r c w (zo;r c ) 1_fe , then Eq. (fT4"]) can be 
simplified as, 
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1-k' 



(15) 



which is the mathematical expression of the Cole-Cole 
equation 0, H[ . 

If k! is large enough and cannot be neglected, we can 
directly take advantage of Eq. (|B1|) in Appendix since 
|iwr c | < 1. Then Eq. (|14[) can be re- written as 



e(ca) 



(1 + iujT c ) 



1-k' 



(16) 



which is the mathematical expression of the Cole- 
Davidson equation 0, 0] ■ 

In Eq. (fTB")) , if we further assume that there exists k" 
and k" <C 1, then 1 — fc" « 1. We then have iwr c w 
i(l — k")uiT c . By using the same approach exploited in 
deriving the Cole-Cole equation, we can get »wr c « i(l — 
k")ujT c w (iu)T c ) 1 ~ k . Substitute this expression into Eq. 
(I1GI) . we have the following relationship: 
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(17) 



which is the mathematic expression of the Havriliak- 
Negami equation 0, Q . 

Case II - the applied electric field is large enough to be 
able to disturb the nematic phase and the term (1 — k) 
is given in the Gibbs free energy. By using the same 
deriving method, we can get the following result. 
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where C = cxp[fc]. This equation is the mathematic 
expression of the Kohlrausch function or the stretched 
exponential function [l(|. Therefore, the Kohlrausch- 
Williams- Watts function (the Fourier transform of the 
Kohlrausch function), which was first used by Williams 
and Watts fll) . can be simplified as 
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(20) 



where 



i-fc- 



Concluding remarks - it is perhaps safe to say now that 
the key factor making the Debye law deviate strongly 
from the relaxation behavior of most dielectric materi- 
als is that the cooperative movement of polar HTSPs 
has been neglected in previous studies. In view of what 
have been derived and discussed above, we believe that 
there are two dynamic processes involved in dielectric re- 
laxation, the first corresponds to relaxation of crystalline 
lattice and the second represents slow relaxation of a par- 
tially ordered phase of polar HTSPs in solid dielectrics. 
The universal dielectric relaxation law does exist. All di- 
electric relaxation relationships including the Dcbye law 
are just the variants of this universal law under different 
situations. 
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Appendix A: Binomial series expansion 

For a function f(x) = (1 + xp, its binomial series 
expansion can be written as 



... + «*- 1 >'-e-" +1 W.., 



Since t £ (0, r c ) and < 1, we can write (1 — k) \Jj- 

k = 1 + (1 - k) (± - l). Obviously, |-|L - 1| < 1. By 
taking advantage of Eq. 
(l-k)(±) + k^(±) 



taking advantage of Eq. (|Blj) in Appendix, we can get 

1-k 

Then the above equation 

can be rc- written as 



Ccxp 



1-k' 



(19) 



here, for simplicity, we only consider the case, in which £ 
is a real number. If \x\ < 1, the above series will converge 
absolutely for any number £ and, thus, can be simplified 
to the following linear equation. 



+ PS 1 + far. 



(Bl) 
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